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Robust pole placement with Moore's algorithm 

Robert Schmid, Amit Pandey and Thang Nguyen 

Abstract 

We consider the classic problem of pole placement by state feedback. We adapt the Moore eigen- 



structure assignment algorithm to obtain a novel parametric form for the pole-placing gain matrix, and 
introduce an unconstrained nonlinear optimization algorithm to obtain a gain matrix that will deliver 



robust pole placement. Numerical experiments indicate the algorithm's performance compares favorably 
against several other notable robust pole placement methods from the literature. 



J. Introduction 

We consider the classic problem of pole placement for LTI systems in state space form 

x(t)=Ax(t)+Bu(t), (1) 

where, for all t £ R, x(t) E MJ 1 is the state, and u(t) G M m is the control input. A and B 
are appropriate dimensional constant matrices. We assume that B has full column rank. We 
let C = {Ai, . . . , A^} be a self-conjugate set of n complex numbers, with associated algebraic 

multiplicities M. = {mi, . . . ,m u j satisfying mi + h m v = n. The problem of exact pole 

placement by state feedback (EPP) is that of finding a real matrix F such that the closed-loop 
matrix A + BF has non-defective eigenvalues in C, i.e F satisfies 



(A + BF)X = XA (2) 

X' 

where A is a n x n diagonal matrix obtained from the eigenvalues of C, including multiplicities, 
and X is a non-singular matrix of closed-loop eigenvectors of unit length. If (A, B) has any 
uncontrollable modes, these are assumed to be included within the set C. The EPP problem has 
been studied for several decades, and the existence of such a matrix yielding diagonal A requires 
the nii to satisfy certain inequalities in terms of the controllability indices of the pair (A, B) [|2l; 
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in particular < m for all m ; e M. is required. In this paper we shall assume (A, B, £, .M ) are 
such that at least one F exists that yields diagonal A. Notable early papers offering algorithms 
to obtaining the required gain matrix F include [3], which gave a method for single-input single 
output (SISO) system, but this was often found to be numerically inaccurate. Varga [0) gave a 
numerically reliable method to obtain F for multiple-input multiple output (MIMO) systems. 

For SISO systems, F is unique, while for MIMO systems it is not, and this naturally invites 
the selection of F that achieves the desired pole placement and also possesses other desirable 
characteristics, such as minimizing the control input amplitude used, and improving numerical 
stability. In order to consider optimal selections for the gain matrix, it is important to have 
a parametric formula for the set of gain matrices that deliver the desired pole placement, 
and numerous such parameterizations have appeared. Bhattacharyya and de Souza A3 gave a 
procedure for obtaining the gain matrix by solving a Sylvester equation in terms of a n x m 
parameter matrix, provided the closed-loop eigenvalues did not coincide with the open loop 
ones. Fahmy and O'Reilly (61, gave a parametric form in terms of the inverses of the matrices 
A — AjJ, which also required the assumption that the closed loop eigenvalues were all distinct 
from the open loop ones. Kautsky et al [7] gave a parametric form involving a QR-factorization 
for B and a Sylvester equation for X; this formulation did not require the closed-loop poles to 
be different from the open-loop poles. 

The classic eigenstructure assignment algorithm of B.C. Moore BH quantified the freedom to 
simultaneously assign both the closed-loop eigenvalues, and also select the associated eigenvec- 
tors. As such it implicitly solved the EPP problem, but it did not explicitly provide a parametric 
formula for the pole-placing matrix, nor did it address any optimal pole placement problem. In 
this paper we adapt Moore's algorithm to obtain a simple parametric formula for the pole-placing 
gain matrix, in terms of an n x m parameter matrix. The method obtains the eigenvector matrix 
X by selecting eigenvectors from the nullspaces of the system matrices, and thus avoids the 
need for coordinate transformations. 

The robust exact pole placement problem (REPP) involves solving the EPP problem and also 
obtaining F that renders the eigenvalues of A + BF as insensitive to perturbations in A, B and 
F as possible. Numerous results [10] have appeared linking the sensitivity of the eigenvalues to 
various measures of the conditioning of X, in terms of the Euclidean and Frobenius norms. This 
classic optimal control problem also has an extensive literature, and typically two approaches 
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have been used to obtain good robust conditioning. 

Perhaps the best-known method for the REPP is that of Kautsky et al 0, which involved 
selecting an initial candidate set of closed-loop eigenvectors and then using a variety of heuris- 
tic methods to make these vectors more orthonormal. This method has been implemented as 
MATLAB®'s place command; this implementation includes a heuristic extension to accom- 
modate complex conjugate pairs in C. This algorithm is also the basis of MATHEMATICA®'s 
KNVD command. The use of the place algorithm has become wide-spread in the control systems 
literature, and introductory texts advocating its use include ifTTI and lfl2l . among many others. 

Since the publication of [7], many alternative methods have been proposed for the REPP. Tits 
and Yang [fT3l revisited the heuristic methods of flT] and offered a range of improvements; the 
algorithms were shown to be globally convergent. Byers and Nash lfl4l . Tarn and Lam [15J and 
Varga [fT6l cast the problem as an unconstrained nonlinear optimization problem, in terms of 
the Frobenius conditioning, to be solved by gradient iterative search methods. [fTTl introduced 
a method for minimizing the 'departure from normality' robustness measure, which considers 
the size of the upper triangular part of the Schur form. Ait Rami et al [fT8| introduced a global 
constrained nonlinear optimal problem in terms of a Sylvester equation and showed that the 
solution could be approximated by a convex linear problem for which the authors gave an 
LMI-based algorithm. 

Various authors have provided surveys comparing the performance of several of these algo- 
rithms. Sima et al |[T9l conducted testing of the algorithms from 0), [0 and [fT3l on collections 
of systems of varying dimensions; they concluded that the method of |[T3l generally gave superior 
Euclidean (2-norm) conditioning and also improved accuracy. [fTTl considered the eleven bench- 
mark systems in the Byers-Nash collection (see Section [IV] for a discussion of this collection), 
and compared the author's proposed methods, based on the Schur form of the open loop systems, 
with those of and [fT3l against a range of robustness measures. The methods of [fTTl generally 
gave inferior results to those of []7| and [fT3l . with respect to the Frobenius conditioning. [fT8l 
tabulated figures results for the Frobenius conditioning performance of methods [Q, lfl3l . Ifl4l 
and [16]. However, the conditioning values were compiled directly from these papers. Since some 
of these methods were introduced into the literature more than two decades ago, and noting that 
computational resources have improved dramatically over this time, using values from original 
publications may unfairly disadvantage the earlier methods, in particular [[141. 
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In this paper we add to this extensive literature in several ways. In Section 2 we introduce our 
parametric form for the pole-placing gain matrix that solves the EPP. The formula is an adaptation 
of the pole placement method of Moore [9|; the novelty here is to use Moore's method to obtain 
a parametric formula for both X, the matrix of eigenvectors and F, the pole-placing gain matrix. 
We further show the parametric form is comprehensive, in that it generates all possible X and F 
that solve ©, for the case where the eigenvalues have multiplicity of at most m. In Section 3 we 
utilize this parametric form to propose an unconstrained optimization problem to seek solutions 
to the REPP, to be solved by gradient search methods. Our approach most closely resembles 
that of |fl4|. but with a different parametric formulation for the pole-placing gain matrix. 

In Section 4 we select five of the most prominent methods for the REPP 0, lfT3l . lfT6ll . 
Ifl4l and [fT8l . and conduct extensive numerical testing to compare their performance against 
our method. The first three of these were chosen as they are widely used in the forms of the 
MATLAB® toolboxes place, robpole and sylvplace respectively. [14] has attracted a large number 
of citations over more than two decades, and [[181 is the most recent publication to offer a novel 
approach for the REPP. All methods were implemented in MATLAB® 20 12a, running on the 
same computing platform. In addition to conditioning, we also compare their accuracy, matrix 
gain and runtime. Finally, Section 5 offers some conclusions as to the relative performance of 
these six methods; our method will be shown to offer some performance advantages over all the 
other methods surveyed. 

II. Pole placement via Moore's algorithm 

We now revisit Moore's method [9] and adapt it to give a simple parametric formula for a 
gain matrix F that solves the pole placement problem, in terms of an arbitrary real parameter 
matrix. We begin with some definitions and notation. For each % G {1, . . . , u}, we define the 
n x (n + m) system matrix 

S(\i) = [A- XJ n B\ (3) 

where I n is the identity matrix of size n. We let Tj be a basis matrix for the nullspace of S(Xi), 
we use Si to denote the dimension of this nullspace, and we denote T =: [T% . . . T u ]. It follows 
that Si = m, unless is an uncontrollable mode of the pair (A,B), in which case we will 
have Si > m. Let M denote any complex matrix partitioned into submatrices M = [M\\ . . . \M V \ 
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such that any complex submatrices occur consecutively in complex conjugate pairs. We define a 
real matrix Re(M) of the same dimension as M thus: if Mj and M i+1 are consecutive complex 
conjugate submatrices of M, then the corresponding submatrices of Re(M) are |(Mj + M i+ i) 
and jj(Mi — M i+1 ). Finally, for any real or complex matrix X of with at least n + m rows, we 
define matrices W(X) and tt(X) by taking the first n and last m rows of X, respectively. 

Proposition 2.1: Let the eigenvalues {A l5 . . . , A^} be ordered so that, for some integer s, the 
first 2s values are complex while the remaining are real, and for all odd i < 2 s we have 
Aj+i = Aj. Let K : = diag(i£i, . . . , K v ), where each Ki is of dimension s« x rrii, and for all odd 
i < 2s, we have K { = K i+1 . Let M(K) be an (n + m) x n complex matrix given by 

M(K) = TK (4) 



and let 



X{K) = 7f(M(K)), 
V(J0 = 7T(Re(M(K))) 
W(K) = n(Re(M(K))) 



(5) 
(6) 
(7) 



For almost every choice of the parameter matrix K, the rank of X is equal to n. The set of all 
m x n gain matrices F satisfying © is parameterised in K as 

F{K) = WiK)V(K)- 1 (8) 

where K is such that rank(X(K)) = n. 

Proof: For any given let M(K) be partitioned according to 

W ... V 
W[ ... 

where each V- and are matrices of dimensions n x m; and m x respectively, such that 

(A-AiW + BWi' = (10) 

Note that, for odd z < 2s, we have that V( = V- +1 are conjugate matrices, as Ki = K i+1 . 
Moreover, since C is symmetric, we also have rrii = m i+ i. Define real matrices 



M(K) 



(9) 



U V i+ V Ui) if«<2sisodd, 
Tj ( V i-i - V i) if « < 2 s is even, 
V' i>2s 



(ID 



DRAFT 



6 



and define Wi similarly. Then matrices X, V and W in ©-([7]) may be written as 

X=[V{ V{ ... vi s \vi s+l V 2 ' s+2 ... VI], V =[V 1 v 2 ... v 2s \v 2s+1 v 2s+2 
waAW=[W 1 W 2 ... W 2s \W 2s+l W 2s+2 ... W v \. Let 



Ri. 



1 

2 



(12) 



I mi 3 Inn 

Then for each odd i < 2s, we have [V{ V( +l ]Ri = [Vi V i+1 ] and \W( W( +1 ]Ri = [Wi W i+1 \. 
Now assume K is such that rank(X(F)) = n; then V(K) is non-singular, and we can obtain 
F in ©. We obtain F[V[ V( +1 ] = [W{ W( +1 ] for odd i £ {1, . . . , 2 s} and FK/ = W/ for all 
i £ {2s + 1, . . . , u}. Hence (flOl) can be written as 



(A + BF) 



V( V' x 



diag(AiJ mi , A i+ iJ mi ),for odd i £ {1, . . . ,2sXl3) 



{A + BF)V( = V l '(X l I m J,for i E{2s + l,...,u}, 



(14) 



Thus we obtain ©. To see that this formula is comprehensive, we let F be any real gain matrix 
satisfying ©. The nonsingular eigenvector matrix X is comprised of column vectors V- of 
dimension n x m ; corresponding to each eigenvalue, such that (fl"3T ) and (fl4T) hold. Applying 
F[V[ Vl +l ] = [WI W( +1 ] for odd % £ {1, . . . , 2 s} and = Wi for alU £ {2s + 1, . . . , u}, we 
obtain V/ and W[ such that (flOl) holds. Thus each column vector of the matrix \V( W[\ T lies 
in the kernel of S(Xi), and we have a coefficient vector Ki such that \V( W{) T = TiKi. The 
complex conjugacy of V( and V- +1 , for each odd i £ {1, . . . , 2 s}, implies the conjugacy of 
and K i+l . Thus we obtain M(K) in © yielding F in ©. 

Finally we let F be arbitrary parameter matrix and consider the rank of X(K). We introduce 
$ = 7f(T) and denote $ l5 . . . , $^ as a basis for im If rank(X(if)) is smaller than n, then 
one column of the matrix [&iK 1;1 . . . § v K v m J\ is linearly dependent of all the remaining ones. 
(Here we have used K^j to denote the j-th column of Ki). For brevity, let us assume this is 
the last column. Then there exist n — 1 coefficients ck 1)1; . . . , a v>mv -\ (not all equal to zero) for 
which 

v— 1 rrii m„-l 
i=l j=l j=l 

has a unique solution in K vmv . As K vmv is an s^-dimensional parameter vector, (fTTb constrains 
^,m„ to lie upon an (sj, — 1) -dimensional hyperplane, which has empty interior. Thus the set of 
parameters K that lead to a loss of rank in X(K) is given by the union of at most n hyperplanes 
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of empty interior. This set therefore has empty interior, and thus also zero Lebesgue measure. 
Thus we see that X(K) and hence V(K) are non-singular for almost all choices of the parameter 
matrix K. ■ 

The above formulation takes its inspiration from the proof of Proposition 1 in O, and hence 
we shall refer to ©-dD as the Moore parametric form for X and F. We note however that 
flU only considered the case of distinct eigenvalues, and did not offer any explicit parametric 
formula for the pole-placing gain matrix. Moreover, it did not show that all matrices X and F 
solving © could be parameterized in the above manner. 

It is interesting to compare this parametric form with that of 0, in which the eigenvectors 
comprising X were obtained from the nullspaces of the matrices U\{A — XJ), where the 
parameter U\ was obtained from the QR-factorization for B = [Uq U\]{Z 0] t , and was also 
required to satisfy U\(AX — XX) = 0. By contrast, the Moore parametric form obtains the 
eigenvectors directly from the nullspaces of the system matrices [A — X-Jn B]. 

III. Robust and minimum gain pole placement 

When A + BF has n distinct eigenvalues, the sensitivity of an eigenvalue Aj of A + BF to 
perturbations in A, B, and F can be represented by the condition number [flOl 

IMI2INH2 



\yf x i 



(16) 



where and xt are the left and right eigenvectors associated with Aj. For the case where A+BF 
is non-defective but has repeated eigenvalues, see ll20l for a definition of the corresponding 
condition numbers. Furthermore, we have [|7] 

Coo := maxQ < K 2 (X) < Kf ro {X) (17) 

i 

where K%(X) = ||X|| 2 ||X _1 1| 2 and Kf ro (X) = ||X||/ ro ||X _1 \\f ro are the condition numbers of 
the matrix of eigenvectors X with respect to the Euclidean and Frobenius norms. Following 
CHI, |[T4|. lfT31 . we propose to address the REPP problem by minimizing the condition number 
of X with respect to Frobenius norm. The objective function to be minimized is 

h{K) = K fro (X(K)) = WXWWfroWX-WWfro (18) 

where the input parameter matrix K is defined as in Proposition 12. II Note it is possible to reduce 
the Frobenius norm of a matrix X by suitably scaling the lengths of its column vectors. When X 
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is the solution to ©, such scaling does not improve the eigenvalue conditioning in (fT6l) . Hence 
we assume that the column vectors of X have been normalised. 

As pointed out in [14J, for efficient computation we can study an alternative objective function 

f 2 {K) = \\X{K)\\) ro +\\X-\K)\\) ro (19) 

because the two objective functions are equivalent. An imported related problem is that of 
minimizing the norm of the gain matrix F. The minimum gain robust exact pole placement 
problem (MGREPP) involves simultaneously minimizing both the conditioning and the matrix 
gain via the weighted objective function 

f 3 (K) = aK fro (X(K)) + (1 - a)\\F(K)\\ fro (20) 

where a is a weighting factor, with < a < 1. Minimizing f 3 involves a gradient search 
employing the first and second order derivatives of Kj ro (X(K)) and ||F(i^)||^ ro ; expressions for 
these were given in [fl]|. 

IV. Performance comparison of robust pole placement methods 

In this section we conduct extensive numerical experiments to compare the performance of our 
method against those of 0, lfT8l . [fT3l . Ifl4ll and lfl6l . To provide a comprehensive contemporary 
survey, we implemented these algorithms on the same modern computer, an Intel® Core™ Quad 
CPU, Model Q9400 at 2.66 GHz with 3326 MB of RAM running Windows™ XP and MATLAB® 
2012a. Implementation of was done with MATLAB®'s place command. For [fT3l and lfT6l . 
we used the robpole and sylvplace MATLAB® toolboxes, kindly provided to us by the authors. 
For [fl4l. ITT8TI and our own method, we wrote MATLAB® toolbox implementations for each. 
The lfl8l algorithm requires an LMI solver; we chose the public-domain cvx toolbox EH . We 
shall refer to these as byersnash, rfbt and span (our own method). The names are derived from 
the names of the respective authors. 

To obtain a fair comparison between these methods, we need to consider the runtime allocated 
to them. The methods of lfT4l . [fT6ll and our proposed method all employ gradient iterative 
searches, so the values they deliver are contingent upon the initial condition (input parameter 
matrix K) used. The sylvplace toolbox randomly generates an initial condition, and thus offers 
different outputs (different F) each time it is run. To obtain repeatable results, we provided 
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the byersnash and span toolboxes with a pre-specified collection of input parameter matrices K 
composed of canonical vectors. The output shown from each of byersnash, sylvplace and span is 
the best result from all the initial conditions searched within the allocated runtime. By contrast 
place, robpole and rfbt all employ a designated starting point, and hence their runtime is simply 
the time taken to execute their method. 

A. Robust conditioning comparison using the Byers and Nash benchmark examples 

Byers and Nash [[T4| gave a collection of eleven benchmark example systems, and many 
authors, including [[T3l . Ifl6l and lfT8l used these examples to compare the performance of their 
pole placement methods. Following this well-established tradition, our first set of comparisons 
employs these well-known examples. The results are given in Table HI We have used Kf ro (X) 
as the performance measure, and we also show the matrix gain used. 

The average runtimes for place, robpole and rfbt for the 11 sample systems were 0.05, 0.095 
and 14.1 seconds, respectively. For byersnash, sylvplace and span we arbitrarily set the runtime 
to be n seconds, where n is the system dimension, leading to average runtimes of 4.5 seconds, 
this being the average of the system dimensions in the collection. 

Ignoring differences in the conditioning of smaller than 1%, we conclude that byersnash and 
span had the best or equal best conditioning in all 1 1 examples, sylvplace and rfbt had the best or 
equal best in 7 cases, while robpole had best or equal best in 5 cases. Finally place gave the best 
or equal best in 4 cases, place and robpole had the shortest runtimes, while rfbt had noticeably 
the longest. We note that the conditioning numbers given here differ significantly from those 
that were published in lfl4l and lfT8l . This may be explained by the fact that these authors did 
not require the columns of X to be of unit length. Since methods and lfT3l normalise the 
columns of X, this is essential for a fair comparison of all six methods. 

B. Robust conditioning comparison with sets of higher-dimensional systems 

To probe more deeply into the performance delivered by these six methods, we need to 
move beyond the low-dimensional examples in the Byers and Nash collection. In Survey 2 we 
generated three sets of 500 sample systems with (A,B), all of state dimension n = 20, and 
with control input dimensions of m = 2, m = 4 and m = 8. The pole positions C were chosen 
to be all distinct, with a mixture of real and complex values. The entries of A, B and C took 
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uniformly distributed values within the interval [—2, 2]. To compare the conditioning, accuracy, 
and matrix gain of each method, we computed, for each system j E {1, . . . ,500} and each 
method * G {place, robpole, byersnash, sylvplace, rfbt, span}, 

• Kf ro (*,j): the Frobenius conditioning of method * for the j-th system; 

• Coo(*, j)' the Cqo conditioning of method * for the j-th system; 

• A(*, j): the accuracy of method * on the j-th system, equal to the largest absolute value 
difference between each eigenvalue of A + BF and the corresponding A, in L. 

• ll-^1l/ro(*} j)' tne Frobenius norm of F from Method * on system j. 

Noting that place is the industry standard for the REPP, we chose to compare all the other 
methods according to their ability to improve upon place, and computed comparative performance 
indices relative to place for each method, and for each performance criterion, as follows: 

500 / .x 

(l-index(*, K/ro )) 500 = TT Kfr .°)*' J) .. (21) 

500 , .s 

(l-mdex(*, Coo )) 500 = ( 22 ) 

500 . , .x 

(l-i„dex(,,A)r = n^T) (23) 



3=1 
500 



:i-i»ta(*,iiF| W r = n nw%±^ < 24 > 



F|| fro (place,j) 



For example, in (|24l) . if mdex(robpole, \\F\\f ro ) = 0.1, then Method robpole gives values of 
||F||j ro that are on average 10% smaller than place. Larger indices imply greater improvement 
on place, and negative indices indicate performance inferior to place. The local gradient search 
methods span, byersnash and sylvplace were each given 20 seconds of runtime per sample 
system; the results shown in Table [IT] represent the best conditioning performance achieved from 
all the initial conditions searched within that time period. For robpole and rfbt, the average 
runtime per sample system were 0.552 and 125 seconds (m = 2), 0.552 and 82.9 seconds 
(m = 4), and 0.552 and 55.2 seconds (m = 8). 

The results show that the best performance for robustness and gain minimisation were given 
by span, byersnash and sylvplace. Both sylvplace and rfbt were less accurate than place, by 
several orders of magnitude in the case of rfbt, which also required substantially longer runtime. 
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While all methods offered improved conditioning with reduced gain over place, this was reduced 
for the larger values of m, which may be attributed to the improved performance of place when 
it has more control inputs to work with. 

C. Weighted gain minimisation and conditioning problem 

Among the methods in our survey, only [|T6l {sylvplace) considered the MGREPP problem 
(|20|) . Our Survey 3 compares the performance of sylvplace and span for the same 500 sample 
systems used in Survey 2, with m = 2, for several different values of the weighting factor a. 
We again gave span and sylvplace 20 seconds of runtime per sample system, and computed 
the performance improvement indices (|2T1 i-d24T) relative to the gain matrix delivered by place; 
again larger figures indicate greater improvement. The results are shown in Table [TTTJ Both 
methods were able to offer significant reductions in gain, at the price of some reduction in the 
robustness measures, relative to the pure robustness problem (a = 1). However span did so with 
far superior accuracy. Considering the impact of different values of the weighting factor, we see 
that for a = 0.1, there was little difference in the conditioning, and only slight improvement in 
the matrix gain. For a — > we observed up considerable reduction in the matrix gain, but this 
eventually comes at the cost of significantly inferior conditioning. These results suggest values 
around a = 0.001 can give a good balance between these two criteria. 

D. Systems with uncontrollable modes 

The EPP problem remains well-posed for systems with uncontrollable modes, provided these 
are included within the set C The methods place, sylvplace, robpole, rfbt all assumed control- 
lability of the system, as part of their problem formulation. In principle this involves no loss of 
generality, since the application of a Householder staircase transformation can decompose any 
system into its controllable and uncontrollable parts. Nonetheless is it is interesting to consider 
the ability of these toolboxes to accommodate uncontrollable modes. In our final survey, we 
obtained 100 systems (A, B), with n = 3 and m = 2, that contained one uncontrollable mode. 
We then chose C to include this mode, plus one pair of complex conjugate modes. We defined 
failure to solve the EPP as being any one of (i) an error was returned upon execution of the 
algorithm, (ii) any of the closed-loop poles differed by more than 5% from their desired location, 
and (iii) the gain of F was undefined or greater than 10 10 . We observed failures as follows: place, 
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sylvplace, robpole and rfbt had 100, 98, 30 and 12 failures, respectively; we concluded these 
toolboxes in their present form cannot reliably solve the EPP in these conditions, byersnash 
and our method span had no failures; we attribute their superior reliability to their usage of 
nullspace methods. Uncontrollable modes increase the column dimension of the corresponding 
nullspace basis matrix; for byersnash and span this is readily accommodated by adjusting the 
row dimension of the parameter matrix. 

V. Conclusion 

We have introduced a parametric formula for the exact pole placement of linear systems via 
state feedback, derived from Moore's classic eigenstructure method. This parametric form was 
used to formulate the robust and minimum gain exact pole placement problem as an unconstrained 
optimization problem, to be solved by gradient iterative methods. 

The method was implemented as a MATLAB® toolbox called span, and its performance 
was compared against several other methods from the classic and recent literature. All methods 
considered gave superior performance to the widely used MATLAB® place command, albeit 
with somewhat longer runtime. When the Frobenius conditioning of the eigenvector matrix is 
used as the robustness measure, the best performance was provided by the our proposed method, 
and also the Byers-Nash method. The results suggest that, in comparison with heuristic methods, 
gradient iterative methods are best able to take advantage of the high levels of computational 
power that are now widely available. They also suggest that methods based on nullspaces of 
appropriate system matrices may offer superior accuracy of pole placement to those adopting 
Sylvester matrix transformations. 

For a given system (A, B, C, M), byersnash and span will in general yield quite different 
gain matrices, offering different performance values, so both methods should be considered for 
optimal performance. While Byers and Nash considered only the robustness, our method is able 
to accommodate a combined robustness and gain minimization approach, enabling the designer 
to obtain significantly reduced gain in exchange for somewhat inferior conditioning. 

The authors would like to thank Andre Tits and Andreas Varga for providing us with copies 
of their robpole and sylvplace toolboxes, and Ben Chen for bringing the classic eigenstructure 
assignment paper by B.C. Moore [HI to our attention. We also thank the anonymous reviewers 
for some constructive suggestions. 

DRAFT 



13 



References 

[1] R. Schmid, T. Nguyen and A. Pandey, Optimal Pole placement with Moore's algorithm, in Proceedings 1st IEEE Australian 

Control Conference (AUCC 2011), Melbourne, Australia, 2011. 
[2] H. H. Rosenbrock, State-Space and Multioariable Theory. New York: Wiley, 1970. 

[3] J. Ackermann, Der Entwurf Linearer Regelungsysteme in Zustandsraum, Regulungstech. Prozess-Datanverarb, vol. 7, 
pp. 297-300, 1972. 

[4] A. Varga, A Schur Method for Pole Assignment, IEEE Transactions on Automatic Control, vol. 26(2), pp. 517-519, 1981. 
[5] S.P Bhattacharyya and E. de Souza, Pole assignment via Sylvesters equation, Systems & Control Letters, vol. 1(4), 
pp. 261-263, 1981. 

[6] M.M. Fahmy and J. O'Reilly, Eigenstructure Assignment in Linear Multivariable Systems-A Parametric Solution, IEEE 

Transactions on Automatic Control, vol. 28, pp. 990-994, 1983. 
[7] J. Kautsky, J. N.K. Nichols and P. Van Dooren, Robust Pole Assignment in Linear State Feedback, International Journal 

of Control, vol. 41, pp. 1129-1155, 1985. 
[8] M. Ait Rami, S.E. Faiz, A. Benzaouia, and F. Tadeo, Robust Exact Pole Placement via an LMI-Based Algorithm, 

Proceedings 44th IEEE Conference on Decision and Control, Seville, Spain, 2005. 
[9] B.C. Moore, On the Flexibility Offered by State Feedback in Multivariable systems Beyond Closed Loop Eigenvalue 

Assignment, IEEE Transactions on Automatic Control, vol. 21(5), pp. 689-692, 1976. 
[10] D. S. Watkins, Fundamentals of Matrix Computations, Wiley, 3rd Edition, 2010. 

[11] G. Franklin, J.D. Powell, and A. Emami-Naeini, Feedback Control of Dynamic Systems. 5th Ed, Prentice Hall, 2006. 
[12] R. Stefani, B. Shahian, C. Savant and G. Hostetter, Design of Feedback Control Systems, 4th Edition, Oxford University 
Press 2002. 

[13] A. L. Tits and Y. Yang, Globaly Convergent Algorithms for Robust Pole Assignment by State Feedback, IEEE Transactions 

on Automatic Control., vol. 41(10), pp. 1432-1452, 1996. 
[14] R. Byers and S. G. Nash, Approaches to robust pole assignment, International Journal of Control, vol. 49, pp. 97-117, 

1989. 

[15] H.K. Tarn and J. Lam, Newton's approach to gain-controlled robust pole placement, IEE Proc. -Control Theory 

Applications. Vol. 144(5), pp. 439-446, 1997. 
[16] A. Varga, Robust Pole Assignment via Sylvester Equation Based State Feedback Parametrization, Proceedings IEEE 

International Symposium on Computer-Aided Control System Design, Anchorage, USA, 2000. 
[17] E. Chu, Pole assignment via the Schur form, Systems & Control Letters vol. 56, pp, 303-314, 2007. 
[18] M. Ait Rami, S.E. Faiz, and A. Benzaouia, Robust Exact Pole Placement via an LMI-Based Algorithm, IEEE Transactions 

on Automatic Control, Vol. 54(2) pp. 394-398, 2009. 
[19] V. Sima, A. Tits and Y. Yang, Computational Experience with Robust Pole Assignment Algorithms, Proceedings IEEE 

Conference on Computer Aided Control Systems Design, Munich, Germany, 2006. 
[20] J. Sun, On worst-case condition numbers of a nondefective multiple eigenvalue, Numerische Mathematik. Vol. 68, pp. 

373-382, 1995. 

[21] M.C Grant and S. Boyd, CVX: Matlab Software for Disciplined Convex Programming, available from http://cvxr.com/. 



DRAFT 



14 



TABLE I 

Survey 1 : REPP with the B yers Nash Examples 



Example 


p/flcel7l 


byersnash\ 14 1 


robpole 1131 




Kf ro (X) 


\\F\\fro 


Kfro(X) 


\\F\\fro 


Kf ro (X) 


\\F\\fro 


i 


6.5641 


1.364 


6.4451 


1.4582 


7.3214 


1.3338 


2 


57.491 


301.37 


50.224 


355.19 


52.972 


224.95 


3 


103.18 


105.06 


46.238 


77.215 


55.987 


49.104 


4 


13.431 


9.899 


13.421 


9.4485 


13.421 


9.4462 


5 


146.18 


4.8496 


142.39 


4.5561 


144.78 


5.4168 


6 


6.0018 


21.5 


5.9633 


23.25 


6.0262 


20.197 


7 


12.375 


233.64 


11.302 


326.35 


12.017 


235.08 


8 


36.986 


15.7600 


6.1824 


28.033 


6.1824 


28.599 


9 


28.682 


2356.5 


23.915 


832.22 


23.937 


823.70 


10 


4.0029 


1.4897 


4.113 


5.2687 


4 


1.5174 


11 


14618 


6692.1 


14510 


6580.8 


14510 


6580.7 


p YMTTinl P 
LiAtllllLUG 


svlvplacei4\ 


rfbtdSi 


span 




K fro{X) 


\\F\\fro 


Kfro(X) 


ll-^ll/f-o 


Kfro(X) 


11*11/™ 


1 


6.5997 


1.4662 


6.5595 


1.5253 


6.4451 


1.4582 


2 


50.042 


327.75 


50.185 


361.01 


50.224 


355.17 


3 


45.741 


72.285 


45.772 


73.582 


46.223 


77.146 


4 


13.421 


9.4465 


13.421 


9.366 


13.421 


9.4432 


5 


141.99 


4.8472 


142.82 


4.3963 


142.39 


4.556 


6 


5.9361 


22.474 


6.4086 


14.771 


5.9622 


23.318 


7 


11.353 


271.17 


12.280 


297.85 


11.301 


271.06 


8 


6.1824 


21.827 


9.381 


39.300 


6.1824 


21.102 


9 


24.23 


903.11 


23.925 


884.84 


23.916 


831.23 


10 


4.113 


1.513 


4 


1.5185 


4 


1.517 


11 


16571 


10716 


14475 


6642 


14510 


6581.3 
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TABLE II 

Survey 2: REPP with higher-dimensional systems 



System Dimension 


Metric 


byersnash \ 14 1 


robpole 1131 


sylvplace\4\ 




span 


n = 20, 
m = 2, 
sys = 500 


Kfra(X)(%) 
Coo(%) 
l|f||/r (%) 

Accuracy (%) 


54.670 
62.047 
23.555 
67.356 


9.8815 
10.620 
1.9292 
26.998 


51.938 
59.759 
22.310 
-1.0082 


41.332 
49.447 
14.337 
-46237 


54.603 
61.983 
23.276 
64.344 


n = 20, 
m = 4, 
sys =500 


Coo(%) 
ll^ll/™(%) 

Accuracy (%) 


37.268 
49.418 
15.677 
45.057 


9.150 
9.8601 
4.3745 
23.760 


36.725 
50.226 
15.524 
-65.586 


31.048 
43.374 
11.163 
-169100 


37.264 
49.400 
15.698 
43.034 


n = 20, 
m = 8, 
sys =500 


K/ro (X)(%) 

Coo(%) 
ll^ll/™(%) 

Accuracy (%) 


15.198 
23.271 
3.7940 
18.525 


7.7702 
10.067 
4.7471 
17.8859 


11.745 
20.848 
3.3979 
-44.635 


12.849 
20.840 
1.7034 
-338240 


15.197 
23.236 
3.7860 
16.225 



TABLE III 

Survey 3: MGREPP with higher dimensional systems (n = 20, m = 2, sys =500) 



Metric 


a = 


= 0.0001 


Q 


= 0.001 


a 


= 0.1 


span 


svlvplacel4\ 


span 


svlvplacet4\ 


span 


svlvplacel4\ 


K fro {X){%) 


-25.578 


23.980 


37.641 


41.906 


53.936 


51.699 


Coc(%) 


-13.540 


33.929 


45.966 


51.379 


61.213 


59.465 


\\F\\fro{%) 


50.319 


38.046 


43.577 


37.740 


27.509 


26.404 


Accuracy (%) 


16.992 


-46.326 


57.833 


-16.025 


65.643 


-1.0463 
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